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Abstract. We present interstellar C n (1334.5323 A) and C II* (1335.7077 A) 
column density measurements along 14 Galactic sight-lines. These sight-lines 
sample a variety of Galactic disk environments and include paths that range 
nearly two orders of magnitude in average hydrogen densities (<n(H)>) along 
the lines of sight. Five of the sight-lines show super-Solar gas phase abundances 
of carbon. Our results show that the excess carbon along these sight-lines may 
result from different mechanisms taking place in the regions associated with 
these stars. 
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1. Introduction 

The interstellar medium (ISM) is a mixture of dust and gas and is responsible for interstel- 
lar extinction of starlight. Extinction is produced by both absorption and scattering and 
is generally represented by a plot of A\/ Av (wavelength dependent extinction/extinction 
at V magnitude) against A -1 (inverse of wavelength). The dominant features in the ex- 
tinction curve are the bump centered at 2175 A and the non- linear Far Ultraviolet (FUV) 
rise (Fitzpatrick & Massa, 1986, 1988; hereafter referred to as FM). The strength of the 
bump and the curvature of the FUV rise varies from one sight-line to another. The bump 
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Figure 1. Profile fits to the AA 1239, 1240 A transitions of Mgn for the fourteen lines of sight. 
Dark lines represent normalized flux and red lines represent the model fits. 



in the extinction curves of the Galactic sight-lines is very strong and must be produced by 
an abundant element in the ISM. Since sulphur is believed to be at most weakly depleted 
in the diffuse ISM (Colangeli et al. 2003), and since nitrogen and oxygen are electron 
acceptors, the carrier of the bump must contain one or more elements from the four mem- 
ber set C, Mg, Si, Fe (Draine 1989). Clayton & Martin (1985) noted that the sequence 
Small Magallanic Cloud (SMC), Large Magallanic Cloud (LMC) and Milkyway (MW) 
appear to have an increasing strength of the 2175 A bump and an increasing C/O ratio, 
suggesting that the carrier of the 2175 A bump feature may be carbonaceous. A majority 
of the proposed models for the 2175 A feature are based on carbon: Polycyclic Aromatic 
Hydrocarbon (PAH; Joblin, Leger, & Martin 1992), graphite (Draine 1989), amorphous 
carbon (Duley & Williams 1981), Quenched Carbonaceous Chondrites (QCCs; Sakata et 
al. 1977) and Hydrogenated Amorphous Carbon (HAC; Mennella & Colangeli 1999). 

Carbon is an important element in the ISM and is the fourth most abundant element 
after hydrogen, helium and oxygen. Carbon is the dominant heat source for the diffuse 
ISM. Carbon grains heat the interstellar gas through photoelectric mechanism and re- 
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combination of electrons with the dust grains (Bakes & Tielens 1994). Interstellar cloud 
cooling is through the rotational and fine structure lines CO and C + (Bakes & Tielens 
1994). In order to fully understand the ISM, it is important to know the distribution of 
carbon among its different phases. 

C II is the dominant form of carbon in neutral interstellar clouds (clouds containing 
Hi and H 2 ) because of its ionization potential (11.26 eV). Despite its importance, there 
are only very few abundance measurements of interstellar carbon (Hobbs, York & Oegerle 
1982; Cardelli et al. 1993; Cardelli et al. 1996; Sofia et al. 1997; Sofia, Fitzpatrick & 
Meyer 1998; Sofia et al. 2004). This may be due to the difficulties associated with the 
measurement of C II column density as the ion produces only very strong and very weak 
absorption features in typical neutral sight-lines. The strong C II transitions at 1036.3367 
A and 1334.5323 A are highly saturated and the column density measurements from 
these transitions are associated with large uncertainties. The transition at 2325.4029 A is 
too weak to be detected for many of the sight-lines. The column density measurements 
mentioned in the above works are all from the weak transition at 2325.4029 A. 

Cn transitions are always associated with excited state (Cli*) transitions and more 
often with 13 Cn and its excited state ( 13 Cn*) transitions. Sofia & Parvathi (2009) 
(hereafter referred to as Paper 1) have developed a method to accurately determine the 
C II and C II* column densities by considering the dominant transitions at 1334.5323 
A and 1335.7077 A together with their nearby transitions at 1334.519 A ( 13 Cn), 
1335.649 A ( 13 Cn*), 1335.6627 A ( 12 Cn*) and 1335.692 A ( 13 Cn*). In Paper 1, 
they have used this method to determine the interstellar C n and C II* column densities 
along six translucent sight-lines (sight-lines having Ay > 1). In the present work, we 
have extended this method to 14 sight-lines observed by the Space Telescope Imaging 
Spectrograph (STIS) aboard the Hubble Space Telescope (HST). 



2. Methodology 

The 14 sight-lines chosen in the present study probe a variety of Galactic disk environ- 
ments and include paths that range nearly two orders of magnitude in average sight-line 



Table 1. List of UVBLUE models used 



Target star 



UVBLUE model 



Target star 



UVBLUE model 



HD 75309 
HD 79186 
HD 91824 
HD 91983 
HD 116852 
HD 122879 
HD 157857 



t28000g40p00k2.flx.bz2 
tl4000g40p00k2.flx.bz2 
t22000g40p00k2.flx.bz2 
t20000g40p00k2.flx.bz2 
tl3000g40p00k2.flx.bz2 
t09000g40p00k2.flx.bz2 
t20000g40p00k2.flx.bz2 



HD 185418 
HD 192639 
HD 198781 
HD 206773 
HD 208440 
HD 210809 
HD 212791 



tl3000g40p00k2.flx.bz2 
tl9000g40p00k2.flx.bz2 
t22000g40p00k2.flx.bz2 
tl5000g40p00k2.flx.bz2 
t25000g40p00k2.flx.bz2 
t30000g40p00k2.flx.bz2 
tl8000g40p00k2.flx.bz2 
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hydrogen density (<n(H)>). Amongst all the sight-line parameters, <n(H)> is the one 
that shows the strongest correlation with the elemental abundances for a given set of 
sight-lines. Therefore the depletion along a given line of sight can be accurately repre- 
sented as a function of <n(H)>. The sight-lines in our sample were earlier studied for 
abundances of Mgn, Pn, Mnii, Mil, Cull and Gen by Cartledge et al. (2006). The 
high resolution STIS observations used in the present study are all taken from the archive 
at the Space Telescope Science Institute (STScI). STIS datasets have the advantage that 
they contain other spectral features required for our analysis (e.g. transitions of Mg n 
at 1239 and 1240 A). In these data, each exposure was made using either of the two 
STIS setups; the 0.2" x 0.2" aperture with the E140H grating centered at 1271 A or the 
same grating exposed through the 0.2" x 0.09" aperture. These data have never been 
used for accurately determining interstellar carbon column densities by the method of 
profile fitting. This is because of the intrinsic difficulties in deriving accurate column 
density values from the damping wings of the C n transition. Inaccurate separation of 
the Lorentzian damping wings from the continuum and lack of any cloud (absorbing com- 
ponent) structure information along the sight-line can lead to large uncertainties in the 
result. 

We have accurately separated the damping wings from the continuum by using ex- 
isting stellar models to fit the stellar continuum in the data. Specifically, we employed 
UVBLUE0, a high-resolution library of synthetic spectra of stars covering the ultraviolet 
wavelength range, developed by the Instituto Nacional de Astrofisica Optica y Electronica 
(INAOE) of Pubela (Mexico), the Brera Astronomical Observatory of Milan (Italy) and 
the Bologna Astronomical Observatory (Italy) . The stellar models covering a large range 
of temperatures, surface gravities and metallicities are computed with the Kurucz (2003) 
set of model atmospheres, by means of the ATLAS9/SYNTHE code (Kurucz 2005). In- 
formation from the simbaeH database is used to identify the target star's spectral type 
and a scan over the parameter space about that value is made to find the appropriate 
stellar model. In addition to that, we have modified the stellar model to account for 
stellar reddening, difference in metallicity from the model, rotational velocity of the star 
and the star's radial velocity shifo We thereby obtained a model which can well repre- 
sent the stellar continuum in the data. As the spectrum of HD 91824 has very complex 
structure for its 1334.5323 A transition, we normalized the continuum by supplementing 
the UVBLUE fit with a fifth order polynomial. Details of the UVBLUE models used for 
different lines of sight are given in Table [T] 

The cloud structure information for the sight-lines were obtained by fitting the Mg n 
transitions at AA 1239, 1240 A. The Mgn profile fits for the fourteen lines of sight are 
shown in Figure [TJ We have constructed the absorption line profile using the parameters 
heliocentric velocities of the clouds (v ) , Doppler broadening parameters of the clouds (b) , 



1 http: / / www.inaoep.mx / ~modelos / uvblue / uvblue. html 
2 http: / /simbad. u-strasbg.fr/simbad/ 

3 Steps followed to modify the UVBLUE model are described in detail in Paper 1 
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The sight-lines marked with * show super-Solar abundance of carbon. 
a logio[N(H)] values are obtained from Cartledge et al. (2006). 
'logio <n(H)> values are obtained from Cartledge et al. (2006). 



(xlO 17 

atoms cm -2 ) (atoms cm" 2 ) 



13.38:* 


1.44e+17 


8.998: 


o.i58:5f 


1.26e+15 


6.808: 


6.898:81 


3.53e+16 


8.708: 


n QnO.001 
< - 8y 0.09 


5.06e+16 


8.688: 


2.618:8! 


9.23e+14 


8.408: 


6.708:81 


5.87e+15 


8.498: 


4.898ig 


7.75e+16 


8.318: 


6.438i| 


8.18e+15 


8.408: 


3.468:81 


3.28e+16 


8.108: 


3.368:11 


2.41e+15 


8.388: 


i2.6?:8 02 


8.93e+16 


8.888: 


8-108:28 


4.31e+15 


8.598: 


6.188:18 


5.76c+16 


8.508: 


3.368:81 


1.45e+17 


8.468: 



V. S. Parvathi et al. 





HD 75309 




JLJV^P' i 



HD 210809 




HD 79 186 



200 



HD 91824 



~ 1.5 

1 1.0 

N 

9 0-5 

p 0.0 



400 2 



HD 185418 



1 


iD 192639 






1 / 


i 










HD 157B571J 





200 
Heliocentric velocity (km s" 1 ) 







u 




[ 1 


[ /HD 212791 





200 
Heliocentric velocity (km s"') 



Figure 2. Profile fits to the Cn and Cn* transitions together with their associated 13 C tran- 
sitions. As the 13 C and 12 C transitions are blended, they are difficult to distinguish from each 
other. Dark lines represent normalized flux and red lines represent model fits. 

column densities of the species in the clouds (N) and the shift in the velocity vector due to 
the instrumental uncertainty (v shift)- The atomic constants of the transition are obtained 
from Morton (2003). In order to find the set of parameter values that produces the best 
fit, we have employed a standard IDL routine 'MPFITFUN0 (Markwardt 2009), which 
performs Levenberg-Marquardt least-squares fit to an IDL function. MPFITFUN fits a 
user-supplied model in the form of an IDL function to a set of user-supplied data. Given 
the data and their uncertainties, MPFITFUN finds the best set of model parameters 
which match the data and returns them in an array. We have the component structure 
information from Mgn transitions. Once again we have made use of MPFITFUN to fit 
the C II and C n* transitions, with an additional model parameter 'N star '. This parameter 
N s tar signifies the fraction of Cn existing in the form of Cn*. We have also used the 
terrestrial value for the isotope ratio of 12 C to 13 C. The Cn and Cn* profile fits for the 
fourteen lines of sight are shown in Figure [5J 



4 http: / /cow. physics. wisc.edu/~craigm/idl/mpfittut. html 
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Figure 3. Gas-phase abundance of carbon on a logarithmic scale plotted against the logarithm 
of average hydrogen density. Associated error bars for the data points can also be seen in the 
figure. The solid straight line with an abundance value of 8.46 corresponds to the Solar system 
abundance of carbon (Lodders 2003) and the dashed lines represent associated lcr error bars. 



3. Results and Discussion 



The C II and C n* column density measurements for the 14 sight-lines in our sample are 
listed in Table [H We have also calculated the gas-phase carbon abundances along the 
sight-lines by adopting Lodders' Solar system value of 12+logio(C/H) = 8.46±0.04 (Lod- 
ders 2003). It is assumed that interstellar elemental abundances can be well substituted 
by Solar abundances. Figure [3] shows gas-phase carbon abundance plotted against loga- 
rithm of average hydrogen density. Five of the sight-lines show super-Solar abundance of 
carbon in gas-phase. The deviation in the positive direction from the Solar value is large 
for the sight-lines toward HD 75309 and HD 206773 whereas it can be treated as marginal 
for the sight-lines toward HD 91824, HD 91983 and HD 208440. HD 75309 is located in 
the Vela supernova remnant (Vela SNR) (Nichols & Slavin 2004). The super-Solar carbon 
abundance along this sight-line can be explained by dust in this evolving SNR being first 
heated and then destroyed by shock wave (Vancura et al. 1994). The dust destruction 
leads to the formation of gas. HD 206773 is a member of Cepheus OB2 association (Pan 
et al. 2005). Based on measurements of CO emission, Patel et al. (1998) suggested that 
the bubble was created likely by stellar winds and supernova explosions. From the 21 cm 
data of Hi, Abraham, Balazs & Kun (2000) proposed that a supernova explosion might 
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have occurred as recently as about 2 Myr ago. This explosion might have expanded the 
pre-existing bubble further. The star in its final stages of evolution ejects its outer layers 
by a supernova explosion or planetary nebula. Carbon is one among those heavy ele- 
ments which are likely to be ejected during a supernova explosion. The ejected elements 
get mixed up with the surrounding ISM thereby creating inhomogeneities in elemental 
abundances. Therefore the different mechanisms taking place in the regions surrounding 
HD 75309 and HD 206773 would have favoured the formation of more gas-phase carbon 
along these sight-lines. 

The sight-line towards HD 79186 is interesting as it shows very large deviation from 
the Solar abundance unlike other sight-lines in our sample. This sight-line shows very 
low gas-phase abundance of carbon with very high error bar. From Figure [5] it is evident 
that the blue wings of the C II line is not continuum fitted very well. Our analysis relies 
on continuum fits based on the UVBLUE stellar spectral library. Although the models 
explore a wide range in parameter space, they cannot be expected to always fit our stellar 
continuua exactly. HD 79186 is one of the cases where the continuum of the star differs 
from the UVBLUE models. When the continuum fits are less accurate, the uncertainties 
in the C II column densities are higher. Because of the stochastic nature of the continuum 
mismatch, it is difficult to quantify this uncertainty in the reported errors. However, 
to be consistent with our method, we avoided using additional low order polynomials 
to accurately fit the stellar continuum of HD 79186. A very low gas-phase abundance 
of carbon along the sight-line towards HD 79186 implies the presence of more carbon 
containing dust along this sight-line. 

It is also observed from Figure [3] that the gas-phase abundance of carbon does not 
follow any strict relation with <n(H)>, unlike gas-phase abundances of other elements 
studied by Cartledge et al. (2006) all of which following a Boltzmann function descrip- 
tion of abundance variation. This is an indication of carbon alone possessing some unique 
characteristics. However, we require a large sample of data in order to completely under- 
stand carbon. 
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